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Abstract 

in 

The Wilson approximate dynamics and the Einstein dynamics are compared for binary 
systems. At the second post-Newtonian approximation, genuine two-body aspects are found 
to differ by up to 114%. In the regime of a formal innermost stable circular orbit (ISCO) 
(yf) ' the both dynamics differ by up to 7%. 

Recently Wilson and Mathews [1] proposed a truncated version of the Einstein field equations 
to treat the coalescence of binary neutron stars in a much simplified but still sufficiently precise 
manner. The main idea of this approximation is to neglect the independent ("true") degrees of 
freedom of the gravitational field, i.e. in particular, the full gravitational radiation content. For 
spherically symmetric processes the proposed scheme is identical with the Einstein equations, 
in non-spherically symmetric dynamical situations, even stationary ones, the proposed scheme 
and the Einstein theory only coincide at the first post-Newtonian level of approximation. It is 
perhaps worth mentioning that, in contrast to the Einstein theory, the proposed scheme should 
allow post-Newtonian series expansions to all orders in (integer) powers of 1/c 2 . 

Most recently, Wilson et al. [2] applied the Wilson scheme to the question of instabilities 
in close neutron star binaries and found the remarkable result that general relativity may cause 
otherwise stable stars to collapse prior to merging. In another recent paper Cook et al. [3] 
tested the Wilson scheme for isolated, rapidly rotating relativistic stars. They found a deviation 
from the Einstein theory of at most 5% which they interpreted as very encouraging for a better 
understanding of binary star evolution. 

In this paper we apply the Wilson scheme to point-like binary systems at the second post- 
Newtonian approximation and calculate the periastron advance as well as the orbital period. For 
circular motion, also the dependence of the angular momentum on the orbital angular frequency 
is given. In addition, we use the Hamiltonian of the Wilson scheme to calculate a formal 
innermost stable circular orbit (ISCO) for binary systems. The obtained results are confronted 
with the corresponding results of the Einstein dynamics. 

In the Einstein theory the periastron advance and the orbital period have been calculated by 
Damour and Schafer [4] starting from a second post-Newtonian Hamilton function. In isotropic 
coordinates and in the center of mass system (Pi = — P2 = P), in reduced varibales (p = P//i, 
r = R/GM) the reduced Hamilton function H = H/fj, reads 

H(T,p,v) = IpS-I-^l-SzV-^KS + zV + ^n-p) 2 ] 
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where i/ = fi/M with // = MiM 2 /M and M M] • M>. Mi and M 2 denote the masses of the 
two bodies. The linear momenta of the bodies are Pi and P2, and R denotes a difference of 
their coordinate position vectors, R = Ri — R2; c is the velocity of light. 

The fractional periastron advance per orbital revolution, k, and the orbital period, P, were 
found to be 
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where E is the total center-of-mass energy (numerical value of H ) and where h is the absolute 
value of the reduced angular momentum J/GMfi. 

The Hamilton function in the Wilson scheme is easily obtained as the two-body special case 
of the n-body matter Hamilton function H mat of Schafer [5] (eq. (3.14)). In reduced form this 
Hamiltonian reads (isotropic coordinates), 
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In the test-body limit, v = 0, the two Hamilton functions If and H ma t are identical as they 
should be on reasons of spherical symmetry of the central body. 

The periastron advance of the dynamics eq. (Q) is easily obtained following the path way of 
Ref. [4]. It comes out in the form 
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As one can see from equations 
following structure 

h - Lh^pN 1 L(k 2 P N + k 2 P N \ 
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where k™ pN denotes the ^-independent terms and k™ pN the z/-dependent terms of k. k and 
kmat are different in the ^-dependent terms of the second post-Newtonian order only. Using the 
Newtonian relation between energy and angular momentum, E = (e 2 — l)/2/i 2 , where e denotes 
the eccentricity of the binary orbit, for the same energy and angular monentum, the fractional 
difference between the two periastron advances at the genuine two-body 2pN level reads 
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The application of this expression to the case of equal-mass binaries (u = 1/4) in circular motion 
orbit (e = 0) gives a fractional difference between k and k ma t at the genuine two-body 2pN level 
of about 63%. 

The orbital period in the Wilson scheme is obtained in the form, again following the route 
of Ref. [4], 
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Analogously to the periastron advance, the two orbital periods have the following structure 

P = P N + h P w + pf) + hp^ + P? N ). (9) 

cr c 1 

Also they are different in the ^-dependent terms of the second post-Newtonian order only. 
For the same energy and angular momentum, the fractional difference between the two orbital 
periods at the genuine two-body 2pN level reads 
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where we have used again the Newtonian relation between energy and angular momentum. 
Compared to Ak in the case of e = and v = 1/4, the expression (10) reaches 6% only. 

For circular orbits the angular frequency, u>, is defined through the expression <3?/-P, where 
<E», the angle advance fore one orbital period, is given by $ = 27r(l + k). Taking into account the 
relation between the energy and the angular momentum for circular orbits, for our two dynamical 
situations, Einstein and Wilson respectively, the relation between angular momentum and orbital 
angular frequency turn out to be 
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The fractional difference between h and h ma t at the genuine two-body 2pN level reads, applying 
the same frequency in both cases, 
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For the case of equal masses, v = 1/4, Ah results in the value of 113.7%. 
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In modifying a procedure of Kidder et al. [6] we calculate now a formal ISCO for the second 
post-Newtonian binary dynamics of the Wilson scheme. The idea of the method by Kidder et al. 
was to add to the equations of motion of a test-mass in Schwarzschild spacetime all ,-dependent 
terms of the second post-Newtonian binary equations of motion. 
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We apply this method to the Hamiltonian (Q) which describes the Wilson approximate 
dynamics at the second post-Newtonian approximation. 

The reduced energy H Q of a test body in Schwarzschild spacetime, in isotropic coordinates, 
reads 

H 0/ , l-l/2c 2 r / 7 1 V 4 p 2 , s 

Augmentation of this expression by the second post-Newtonian ^-dependent terms, H u , of the 
Hmiltonian (|3|) yields the following so-called hybrid approximation 

H* = H + H U . (15) 

For circular orbits this expression can be expressed as a function of r and the reduced angular 
momentum h (n • p = 0, p 2 = h 2 /r 2 ). The radius of the ISCO is than obtained by the aid of 
the equations 

dH*(r,h) = Q = 8 2 H4r, h) 
dr dr 2 
We have solved these equations numerically. For two equal masses, v = 1/4, we obtained for 
the ISCO 

r mat = 6.82 GM/c 2 . (17) 

Starting from the second post-Newtonian two-body Hamiltonian (H) the ISCO has been calcu- 
lated by Schafer and Wex [7]. They used the same method and they obtained r = 7.34GM/c 2 . 
The result ( |l7| ) obtained from the Hamiltonian (Q) differs by 7.1% from this value for r. 

For an easy comparison of the two different dynamical situations the numerical values of the 
energy, H, and the angular momentum, h, both coordinate independent quantities, are given in 
Table 1. 



Tab. 1 : ISCO radii r (unit GM/c 2 ) in isotropic coordinates and related angular momentum 
J (unit vGM 2 /c) and energy E (unit uMc 2 ) obtained (i), for a test body in Schwarzschild 
spacetime, TBSD, (ii), for the two-body hybrid 2pN Wilson approximate dynamics, HWAD, 
and (iii), for the two-body hybrid 2pN Einstein dynamics, HED, 





TBSD 


HWAD 


HED 


r 


4,95 


6,82 


7,34 


J 


3,46 


3,59 


3,71 


E 


-0,057 


-0,050 


-0,047 



The results obtained for the periastron advances, eqs. (2) and (5), and the orbital periods, 
eqs. (3) and (8), are valid in dynamical regimes where post-Newtonian appproximations apply. 
In the test-body limit as well as in the first post-Newtonian approximation the periastron 
advances and the orbital periods coincide. They only differ in the genuine binary 2pN parts, 
this means in the ^-dependent terms of 2pN order. The fractional difference of the 2pN in- 
dependent terms for h takes the remarkably large value of about 114% in the case of equal-mass 
binaries in circular motion with the same orbital frequency. 

On the other side, the estimated value for the ISCO of the hybrid 2pN Wilson scheme, eq. 
(14), differs by about 7% from the previously calculated value for the 2pN binary dynamics, 
d), by Schafer and Wex [7]. In Tab. 1 we have summarized the numerical values for the 
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ISCO radius, the energy and the angular momentum for the different models. Surely, near the 
last stable circular orbit post-Newtonian approximations loose their meaning. Nevertheless, it 
seems reasonable to conclude that in the Wilson scheme the binary system, near the ISCO, 
is stronger bounded (see Tab. 1). The conclusion is supported by the fact the in regimes 
where the 2pN approximation applies the stronger binding holds, see Fig. 1 : Given the same 
orbital angular frequency, the Wilson approximate dynamics has smaller angular momentum, 
i.e. smaller moment of inertia. Another support results from the energy expression as function 
of u, most easily obtained from the relation dE = codh, also showing stronger binding. 
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Fig. 1 : The value of the genuine two-body 2pN term in h is plotted as function of u>, the orbital 
angular frequency for circular motion, for the both models. The upper curve belongs to the 2pN 
Einstein dynamics whereas the lower curve belongs to the 2pN Wilson approximate dynamics. 



